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ABSTRACT 

The spectrum of supersymmetric domain wall solitons of the Wess-Zumino 
model is known to be discontinuous across a curve (of marginal stability) in the 
moduli space of quartic superpotentials. Here we show how this phenomenon 
can be understood from the behaviour of the long-range inter-soliton force, 
which we compute by a method due to Manton. 
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1 Introduction 


It is often stated that solitons of supersymmetric field theories that preserve the same 
fraction of supersymmetry will exert no force on each other. When true, this implies 
the existence of static supersymmetric multi-soliton solutions that can be interpreted as 
individual solitons in static marginal equilibrium due to a cancellation of attractive and 
repulsive forces. However, static supersymmetric multi-soliton solutions may not exist. 
This is typically the case for domain walls; i.e. solitons of (l+l)-dimensional field theories. 
Although it is possible to find (l+l)-dimensional models that admit multi-soliton solutions 
of the type described above Q, these seem to be the exception rather than the rule. In 
general, multi-soliton solutions of (l-t-l)-dimensional field theories are time-dependent, 
and this implies the existence of a force between the constituent solitons. 

Consider the case of a field theory for a real scalar field <f> with a potential having 
three isolated degenerate minima, at (f> a < <t>b < 4>c- There exist static soliton solutions 
interpolating between the adjacent minima, and these are supersymmetric solutions of 
the supersymmetric version of this model, but there is no supersymmetric soliton inter¬ 
polating between the non-adjacent minima. This is because all supersymmetric solutions 
correspond to flows of a first-order equation and there is no flow connecting cf) a with <f> c . 
There must be some solutions that interpolate between the non-adjacent minima, but 
they are necessarily time-dependent. These solutions represent an afr-soliton followed by 
a &c-soliton moving under the influence of the force between them, at least for large sep¬ 
aration. The leading-order term in an asymptotic expansion of the long-range force must 
be repulsive because an attractive force would imply the existence of a bound state of 
an afr-soliton with a 6c-soliton and hence the existence of an ac-soliton but, as we have 
just explained, there is no such soliton. We shall confirm the repulsive nature of the 
long-range force, as a special case of a more general result, by adapting a method intro¬ 
duced by Manton || to compute the long-range attractive force between a soliton and its 
anti-soliton. 

The situation for multi-component scalar field theories is much more complicated. 
Here we concentrate on domain walls of the bosonic sector of the Wess-Zumino model 
for a single complex scalar superfield. On reduction to (1+1) dimensions this becomes a 
model for a single complex scalar field z(t, s ) with Lagrangian density 


£ = l [\z\ 2 - \zf - \dW(z)\ 2 ] , 
where the ‘superpotential’ W ( 2 ) is a holomorphic function, and 




dW = 


dW 

dz 


( 1 ) 

( 2 ) 


Critical points of W are degenerate global minima of the potential and solitons are minimal 
energy configurations that interpolate between them. In the supersymmetric context, the 
critical points of W are supersymmetric vacua and the solitons interpolating between 
them preserve 1/2 of the supersymmetry Ill- 

Polynomials provide a simple class of superpotentials. A polynomial superpotential 
of order n has n — 1 critical points, and hence n — 1 vacua, so in order for such a model 
to admit multi-soliton configurations we need n > 4. Here we shall concentrate on the 
simplest case of a quartic superpotential, which may (without loss of generality) be put 
in the form 

W(z) =z 4 - \pz 3 - 2z 2 + 4 pz , (3) 

O 

where p is a complex constant that parametrizes the space of physically-distinct quartic 
superpotentials Q. Provided that p ^ ±1, there are then three (degenerate) supersym- 
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Figure 1: The curve of marginal stability separating the ‘2-soliton’ and ‘3-soliton’ regions 
in the /(-plane, where fi = /ti + i\ii- 


metric vacua, at 


(1) : z = —1, (2): z = l, (3) : z =» , (4) 

and hence, potentially, three soliton solutions interpolating between them. However, 
whether all three types of soliton actually exist depends on the value of /r. For example, 
if /Li = 0 then soliton solutions z(s) are real (because all three vacua lie on the real axis in 
the 2 -plane) and can connect only adjacent vacua; there are therefore only two types of 
soliton. On the other hand, the choice 

/t = iV 3 (5) 

yields a .^-symmetric model for which the existence of all three solitons is guaranteed by 
symmetry given any one of them0. 

The /(-dependence of the soliton spectrum of the WZ model with superpotential (||) 
was analysed in Q, where it was shown that the ‘2-soliton’ to ‘3-soliton’ frontiers in the 
/(-plane are the two branches of a curve A(hi,H 2 ) = 0, where we have set n = /i i + i\in 
and 

A = 3 — 6/4 — 6/4 — a 4 + 2/if/Li| + 3/4 • (6) 

This curve is shown in Fig. 1; when it is crossed from inside a ‘3-soliton’ region, one of 
the three solitons disappears from the spectrun]j|. The curve A = 0 is therefore a simple 
(and, apparently, the earliest) example of what is now called a curve of marginal stability. 
The reason for the discontinuity in the soliton spectrum across this curve was explained in 
Q: consider what happens as we start from the ^-symmetric case (|s|) and proceed down 
the imaginary /r axis to the origin. As /(2 decreases, the soliton trajectory in the 2 -plane 
that starts at vacuum 1 (z = —1) and ends at vacuum 2 (z = 1) passes increasingly close 
to vacuum 3 (2 = /t). In other words, the 12-soliton (connecting vacua 1 and 2) looks 
increasingly like a loose bound state of the 13-soliton and 32-soliton. This suggests the 
following picture: near the curve of marginal stability, one of the three solitons can be 
viewed as a bound state of the other two in which the constituents are held at a distance 
that goes to infinity on the curve of marginal stability, thus causing the bound state 
soliton to disappear from the spectrum. 

3 This case yields an integrable model in (1+1) dimensions H. 

4 This phenomenom was rediscovered by other methods in jf j, where what we call a WZ model is called 
a ‘Landau-Ginzburg’ model. 
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The main aim of this paper is to confirm this picture for the WZ model described 
above by a determination of the asymptotic intersoliton force using Manton’s method j||. 
For real /i all solitons are real and hence solutions of a truncated theory involving only a 
single real scalar field. As explained above, the two solitons of this theory must repel each 
other, at least asymptotically. As we move the z = fi vacuum away from the real axis 
towards the curve of marginal stability we find that this repulsive asymptotic force goes 
to zero, changing sign as we cross the curve. So the leading-order force is attractive on the 
‘3-soliton’ side of the curve, as we might expect. We shall argue, albeit less directly, that 
the next-to-leading order force is always repulsive. This implies the existence of a bound 
state near the curve of marginal stability (on the ‘3-soliton’ side) with a separation of the 
constituent solitons that diverges on the curve. This bound state therefore disappears 
from the spectrum as the curve is crossed, in agreement with the results of Q. 

A similar result was obtained recently, by different methods, for a different scalar field 
theory Jr] . We should also note that similar results have been obtained, again by different 
methods, for monopoles and dyons of (3+l)-dimensional supersymmetric gauge theories 
1 - 110 . 

We shall begin with a brief review of the solitons of the WZ model with quartic 
superpotential. Although they are not known explicitly we will show that the soliton 
trajectories in field space can be found exactly ; this allows some qualitative results of [Q 
to be made quantitative. Next, we show how Manton’s computation can be generalized 
to yield the asymptotic long-range force between solitons of multi-component scalar field 
theories. We then apply this result to the solitons of the WZ model with quartic super¬ 
potential. In particular, we show that the leading-order force vanishes on the curve of 
marginal stability. We conclude with a discussion of the next-to-leading order and the 
implications for soliton bound states. 

2 WZ solitons for quartic superpotential 

Let Zi (z = 1, 2, 3) be the three critical points of the quartic superpotential (j|). Given the 
existence of a soliton connecting vacua i and j, its topological charge is 


Tij = 2 [W( Zj ) - W{ Zi )) . (7) 

These three charges form the three sides of a triangle in the complex W-plane. The 
triangle inequalty 

\Tij\ + \T jk \ > \T ik \ (8) 

ensures that any soliton in the spectrum of supersymmetric states must remain in it as fi is 
varied, unless the triangle degenerates. If this happens, the inequality is saturated and the 
minimum energy configuration with charge need not be a one-soliton configuration. 
In fact, since solitons correspond to straight lines in the W-plane[] only two of the three 
possible solitons can exist whenever the three points W(zi) in the W-plane are colinear; 
this occurs when fi is real and when fi lies on the curve A = 0, with A given by (Q). When 
fj, is real, the long-range intersoliton force is repulsive (for reasons explained above) and, 
by continuity, it remains repulsive for fi in some neighbourhood of the real axis. Thus, in 
this region of moduli space only two of the three possible soliton states exist and there is 
therefore no discontinuity in the soliton spectrum across the real fi axis. This agrees with 
the analysis of Q but it was also shown there, by a study of the qualitative behaviour of 
soliton trajectories in the z-plane, that there is a discontinuity across the curve A = 0. 
We shall now confirm this by finding the exact soliton trajectories. 

5 A recent demonstration of this may be found in fn|. 
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Supersymmetric solitons solve the first order ‘BPS’ equation 

z' = e~ ia dW(z), (9) 

where a = arg(T). Writing z(x) = u{x) + iv{x), this yields the pair of coupled differential 
equations 

u' = f(u,v), v' = g(u,v) (10) 

where, for the quartic superpotential (pj), 

f{u,v) = 4 cos a [u 3 — 3uv 2 — u — gi(u 2 — v 2 — 1) + 2uv/i2\ 

— 4 sin a [u 3 — 3 u 2 v + v + 2uvg\ + ( u 2 — v 2 — 1 )^ 2 ] 

g(u,v) = 4cosa [u 3 — 3u 2 v + v + 2uvg\ + (w 2 — v 2 — l)/i 2 ] 

+ 4 sin a [ u 3 — 3 uv 2 — u — {v 2 — u 2 — l)g.i + 2uv)jl2] • (11) 

These equations determine u(x) and v(x). When /i is not real we cannot find explicit 
solutions, but we can find the soliton trajectories in the z-plane. To this end, we observe 
that & implies that fdv — gdu = 0. However, 

fdv — gdu = dh (12) 


where 
h(u,v ) 


4 cos a 


— sma 


(u 3 — uv 2 — u) v — niv ( u 2 — -v 2 — 1 j + n 2 U ( v 2 — T u 2 + 1 


2v — 2u +v +u — 6 u v + 4/ii u ( v — —u + 1 


+ 4^ 2 f ( u - -v - 1 


(13) 


Soliton trajectories are therefore curves of constant h in the z-plane that pass through two 
of the critical points. All other curves of constant h are ‘BPS-flows’ in that they solve the 
first-order equations ©, but not with the boundary conditions needed for finite energy. 

Consider the 12-soliton interpolating between the vacua at z = — 1 and z = 1. The 
curve of constant h that passes through these points has 



. 1*2 
tan d \2 — — ? 

1*1 


(14) 


The latter relation implies that argT = actan(/i2/Mi)> which is consistent with the fact 
that T oc g, for the 12-soliton. The 12-soliton trajectory is therefore 


0 


M2 


( 1 _ u 2) 2 + i; 2 ( 2 + , 


+ 4|/^| 2 n 




4 fiiuv ( u 2 — v 2 — l) 


(15) 


In the limit as g —> 00 this trajectory coincides with the real axis, but for finite g it is 
energetically favourable for the 12-soliton to partially roll down the potential towards this 
third vacuum. A 12-soliton trajectory that passes very close to z = n can be viewed as a 
bound state of the 13 and 32 solitons. In the limiting case in which the trajectory passes 
through z = fi it must either end or begin there, so we then have an infinitely separated 
13 and 32 soliton but no 12 soliton. Demanding that @ pass through z = fi yields 
precisely the same condition as found in Q by requiring colinear topological charges; that 
is, either g, is real or it lies on the curve A = 0. 
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Figure 2: Plot showing the soliton £12 disappearing from the spectrum as the 2 = /i 
vacuum crosses the curve of marginal stability in the 2 -plane. All plots have /ii = 0 
and fi 2 is respectively 1, 0.7 and 0.65. The critical value for which A = 0 if /ri = 0 is 

A*2 = ±V2Vi -3 « ±0.681. 


The other two soliton trajectories can be found similarly. In each case the soliton 
trajectory coincides with the union of the other two when ^2 A = 0 , although it disappears 
from the spectrum only on crossing the curve A = 0. It should be noted that the relevant 
segments of the curve A = 0 differ in all three cases. The above case, in which the 
12-soliton is the one that disappears from the spectrum as we go from a ‘3-soliton’ to a 
‘ 2 -soliton’ region, corresponds to the segments with fi\ < 1 , whereas the other two cases 
correspond to segments with n\> 1. We skip the details but note for future use that 

tan oi 3 = ——7 — ) tan a 2 3 =-7-7— > ( 16 ) 

fii ± A + /ii ± A_ 

where 

A± = ±—--^=-7 . ( 17 ) 

4 [3/xi ± nml ~ Mi ± 2] 

It follows that oi 3 = 032 when A = 0, as expected. 


3 The leading-order long-range force 

We now plan to obtain a formula for the asymptotic force between two solitons of a (1+1)- 
dimensional field theory for a multi-component real scalar field (f>(t , s ) with Lagrangian 
density 

£=^| 2 -^f-V. (18) 

This of course includes (jl) as a special case of a two-component real scalar field Lagrangian 
density. We assume that V has multiple degenerate vacua at </> = </>;(* = 1, 2,...), and that 
some vacua are connected by static soliton solutions. Let <piji s ) b e a soliton connecting 
the *th vacuum to the jth vacuum. Near the ?’th vacuum we have 

V ~ (19) 

The asymptotic form of the ij-soliton solution near the _)th vacuum (as s —> 00 ) is therefore 

Zij(s) - Zj ~ -mJ 1 t ij e~ rn:iS , (20) 
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where t,j :] is tangent to the zj-soliton trajectory in field space as it approaches the jth 
vacuum. Similarly, the asymptotic form of the zj-soliton near the zth vacuum (as s —> — oo) 
is 

Zij(s) - Zi ^ m~ 1 t ji e miS . (21) 

Now consider a static configuration of an zfc-soliton at the origin, s = 0, separated by 
a distance L from a kj- soliton at s = L. A configuration that fits this description is 


<f>{s) = <f>ik(s) + </>kj(s ~ L) - (f> k . ( 22 ) 

We should not expect this to be an exact solution of the held equations but it will be an 
approximate solution for large L provided that there is a b in the range 0 -C & A such 
that both individual static soliton solutions <j>ik{s) and (p k j (s) are nearly equal to their 
common vacuum value (j> k near s = b. We shall write (|2^) as 

4>{s) = <t>ik(s) + ip(s) (23) 

where 

<f(s) = <t>kj{s - L) - (j) k (24) 

and assume that ip(s) is a small perturbation to the static solution 4>i k (s) in the region 
s <b. Since, for s ~ b, 

<j>ik(s) - <t>k ~ ~m~ l ; 1 ti k e~ rnkS , <p(s) ~ m]; 1 tjke rnk( - s - L ' > , (25) 

this means that we must have 2b <C L , so we are now assuming that 


<f> ■ cf> ds . 


0< 26 < L. 

The force exerted on the ik-soliton by the kj-soliton can now be found as follows 
the momentum of the ik-soliton is approximately given by^j 

— f 

The force on it is therefore 

F = P = 

which follows on use of the held equation 


-\w\ 2 + V 


1 b 


dV 
d 4> 


Using the ansatz (23) and properties of <j>ik , notably 

I<(4| 2 = 2 V{<tH k ), 

we hnd that 


(26) 
s 0: 

(27) 

(28) 

(29) 

(30) 

(31) 


F =~ [4>'ik ■ ¥>' - $lk ■ vl-oo + °W 2 ) ■ 

This may be evaluated using the asymptotic forms (^5|) of <pi k and ip. The result is 

F = -2 t lk ■ t jk e~ mkL . (32) 


This generalizes the result of Manton to multi-component scalar held theories. For a 
single-component theory the tangent vectors are necessarily either parallel or antiparallel, 
so the force is either attractive or repulsive, never zero. For example, for a soliton- 
antisoliton pair we have i = j and hence the attractive force 0 . 

F = -2 t 2 ik e~ mkL . (33) 

6 The sign, which agrees with is chosen such that a negative force corresponds to an attractive 
intersoliton potential. 
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4 Application to the WZ model 

The above results are applicable to the WZ model because this has degenerate vacua at 
critical points z = Zi of the superpotential W, with 

m? = \W'(z i ) | 2 . (34) 

We shall consider again the case of a quartic superpotential with critical points at z = 
— 1,1, fi (with fi = fix + i/j , 2 ), and choose 


Mi < 1 > 


(35) 


so that the 12-soliton is the one expected to appear as a bound state (of a 13-soliton and 
32-soliton) near the curve of marginal stability. Consider a field configuration in which 
a 13-soliton is at s = 0 and a 32-soliton at s = L, for large L. In the region between 
these solitons, and far from both, we may linearize the first-order equations © about 
the z = h vacuum to get 

(36) 


u- hi 
v - H2 


where the matrix M is 


= 4 


Mi - M 2 - 1 2/.ii/i 2 

2hiH2 Mi “ M 2 ~ ! 


cos a 
— sin a 


sma 
cos a 


This matrix has eigenvalues ± 777 - 3 , where 

m 3 = 4y/l + |/i| 4 . 

The corresponding eigenvectors are 

1 


e±{a) = 


\f2 [1 =F cos(0 + a)] 2 


- (sin(a + 9), — cos(a + 6) ± 1) 


where 


tan 9 = 


2hiH2 


(37) 


(38) 


(39) 


(40) 


1 + m! - Mi ’ 

These eigenvectors are tangents to the separatrix flows at z = H- Note that they are 
orthogonal. 

The generic solution of (36) has both e~ m3S and e m3S terms, but the asymptotic 
13-soliton solution 213 has only the e~ m3S exponential term. Thus, we should take 


ti 3 oc e_(aq 3 ). 


(41) 


Similarly, the asymptotic 32-soliton solution, near the 3-vacuum, has only the e m3S factor, 
so we should take 

t 23 oce+(a 32 ). (42) 

The sign of the constant of proportionality should be the same in both cases, so this yields 
the force 


F oc — sin 


(oq 3 — a 32 ) 


Sill 


(oq 3 + a 32 ) + 9 


Sill 


(ai 3 — a 32 ) 


(43) 


for some positive constant of proportionality. 
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Consider first the case for which /Z 2 = 0, so that both soliton solutions correspond 
to BPS-flows along the real axis in the 2 -plane from 2 = —1 to 2=1 (since /zf < 1 by 
hypothesis). From ([h]) we then see that 

OP3 = 0,7T . (44) 


Similarly for CK 32 but 032 = 013 + 7 r. To see why, note that a BPS flow just above the real 
axis will approach the vacuum at 2 = /i anti-parallel to the real axis (e_) and then leave 
it parallel to the imaginary axis (—e+). To arrange for this flow to reach 2 = 1 (or close 
to it) we must rotate —e+ back to e_. This is achieved by a shift of the angle by n since 


e T (a + 7r) = —e±(a). 


(45) 


Thus, in this case, 

£*32 = ai 3 + 7T, (46) 

and this leads to F oc (1 — sin 0). Since sin0 = 0 for real /z we thus deduce that F oc 1 with 
positive constant of proportionality; that is, F > 0, so the force is repulsive, as expected. 

We now turn to those cases for which /z is near the curve of marginal stability. In this 
case we see from © that 


(ai 3 — 023) = 


-A*2 A 


4m?( 1 — An) 3 


+ O (A 2 ) 


(47) 


It follows, after some calculation, that 


F oc 


Ml (1 + l^l 2 ) (l + /z!-/z?) A 


(i - V^a - imI 2 ) 2 +m la + w? (i+ti- Mi)" 


Since /z 2 < 1, all factors other than A are positive. We thus have 

F(L) ~ c 2 [A + O (A 2 )] e" m3i , 


+ O (A 2 ) . (48) 


(49) 


for some non-zero constant c. The leading-order asymptotic force is therefore repulsive 
when A > 0, which corresponds to a point inside the ‘2-soliton’ region. It is attractive 
when A < 0, which corresponds to a point inside a ‘3-soliton’ region. On the curve of 
marginal stability the leading-order asymptotic force is zero. 


5 Soliton bound states 

So far, we have analysed the behaviour of solitons in the WZ model with quartic super¬ 
potential near the curve of marginal stability in the moduli space of such superpotentials. 
On one side of this curve there are only two types of soliton and the long-range asymptotic 
force between these two solitons is repulsive. We have shown, however, that this repulsive 
asymptotic force vanishes on the curve of marginal stability, and becomes attractive after 
it is crossed. 

An interesting question that this result raises is whether the long-range intersoliton 
force vanishes on the curve of marginal stability only to leading order in an asymptotic 
expansion (in powers of e -7 ™ 3 ^) or to all orders (assuming the existence of this expansion). 
If the force were to vanish exactly on the curve of marginal stability then we would expect 
to be able to find a one-parameter family of 12-soliton solutions corresponding to a 13 
and 32 soliton at arbitrary separations. However, there are no such solutions because the 
12 -soliton trajectory necessarily coincides, on (the appropriate segment of) the curve of 
marginal stability, with the union of the 13 and 32-soliton trajectories. This fact suggests 





that it is only the leading-order asymptotic force that vanishes on the curve of marginal 
stability, and that the next-to-leading order term is non-zero. Given that it is non-zero 
it must be repulsive because an attractive force would imply a bound state ‘third soliton’ 
on the ‘2-soliton’ side of the curve of marginal stability. Near A = 0 we thus expect an 
asymptotic expansion for the intersoliton force of the form 

F{L) ~ c 2 Ae~ m3L + y 2 e~ 2m3L + ... (50) 

for some non-zero constant 7 . When A < 0 this force vanishes for 

e-^ L « ( c 2 / 7 2 ) |A|. (51) 

Since A is small the neglect of the higher-order terms in (^0|) is justified, and there is 
thus a minimum of the inter-soliton potential at L ~ ( 1 / 7713 ) log(l/|A|). This diverges on 
the curve of marginal stability, as we know must happen from the analysis of the soliton 
trajectories, and this explains the discontinuity of the spectrum on crossing this curve. 
As mentioned earlier, this result agrees qualitatively with a similar result obtained for a 
different scalar field model by different methods jfj, as well as with results for other types 
of supersymmetric soliton. 
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